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i 

The  ras  dyraraloe  re  ulting  from  the  release  of  an  isothermal 

sphere  of  gas  initially  at  rest  and  at  a  high  pressure  are 

described  by  the  numerical  solution  of  the  differential  equations 

which  represent  an  ideal  gas  in  radial  motion.  The  resulting 

outward  shock  and  the  inward  rarefactions  end  shocks  are  described 

in  sore  detail  and  are  compared  with  the  gas  dynamics  of  the 

point- source  solution,  The  existence  of  tho  inward  moving  shocks. 

2«J< 

-  hae  been  discussed  In  several  places  ^  and  In  one  lnetance  the 

5  V 

early  ohase  has-  been  numerically  calculated,  Thie  pteenomsren  is 
observed  in  the  present  calculations  to  occur  repeatedly,  sub  Joe  ting 
the  origin  to  several  compressive  wav«3  which  are  reflected  there 
and  which  move  out  to  overtake  the  main  shock, 

*H.  L.  Brode,  Phys.  Her.  9$,  656  (A)  (195*). 

^F,  Wecken,  Z,  angew,  Math,  Mach,  ^  270  (1950). 

3J.  A,  HcFadden,  J,  Appl.  Phys.  Q  1269  (1952). 

,4H.  Schardln,  Communications  on  Pure  and  Applied  Mathematics 
(Inst,  of  Mathematical  Sciences,  Mew  York  University,  (195*),  Vol,  VII, 
p.  223. 

5T.  S.  V.'mlton,  Phys.  P«r.  87,  910  (A)  (1952). 
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THB  BLAST  FROM  A  SPHXRB  OP  HZ  OH  FRS8URB  OAS# 

I  hart  UMd  a  method  prcpoaad  in  1950  by  ron  Baumann  and  Ricbtmyer1 
far  numerically  integrating  tha  differential  equations  of  state  and  motion 
of  an  ideal  gas.  TVs  astbod  introduces  an  artificially  large  viscosity  for 
tbs  purpose  of  avoiding  shock  discontinuities  in  tbs  numerics. 

Z  described  tbe  blast  wave  resulting  from  a  point  source  at  tbs  last 
Washington  meeting  •  These  results  are  also  to  be  published  in  tbs  Journal 
of  Applied  Physics.  In  this  paper  I  shall  not  mention  the  details  of  the 
numerical  methods  but  shall  describe  briefly  some  interesting  features  of 
the  blast  wave  resulting  from  a  finite  source,  in  particular  an  initially 
static  and  isothermal  high-pressure  sphere. 


KjUATIOWI  OP  MDTZOM 


The  equations  of  notion  in  some  units  and  in  Lagrangean  form  may  be 

expressed  as  these 


&P 

3; 


.  ^  x  *  *  <^-1) « 


In  which  X  is  a  reduced  radial  distance, 


(energy) 

x  is  the  mass  or  Lagrangean 


Presented  at  the  American  Physical  Society  Nesting,  January  27-29*  1955. 
von  Seumann,  J.  and  Rlchtmysr,  R.  D. ,  J.  Appl.  Phys.  21,  232  (1950) 

p 

Brods,  F.  L. ,  numerical  Solution  of  a  Spherical  Blast  Wave.  Phys. Rev.  658(A) 

1 

Brods,  H.  L. ,  Journal  Appl.  Phys.  (to  be  published) 
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coordinate,  t  it  the  time  coordinate  and  p,  P,  u,  q,  are  respectively  the 
density,  pressure,  particle  velocity,  and  artificial  viscosity  in  units  of 
some  at  >ient  values.  For  this  probleo  gamma  was  taken  to  be  that  of  noraal 
air  (1.4). 

The  particular  form  used  for  the  artificial  viscosity  is  the  followings 


•  -|fj) 

M  »  Ifuraber  of  tones  spread  (2) 

X  -  r/c 


.3-  Etot  4*  f 

-  f —  r 

0  O  Jr 


^int  ♦  r> A* 


1)110  fonr  has  the  property  that  it  la  aero  except  for  compressions,  l.e. , 
shocks,  where  It  spreads  the  ehock  over  a  constant  nusber  of  space  tones  (M). 
The  reduced  radius  (x)  la  shown  to  be  dimensionless  and  related  to  the  total 
blast  energy  (Bt0^)  and  the  ambient  or  pre-shock  pressure  (PQ). 


Dgmogg  BQUATI0B8 

We  approximate  the  differential  equations  by  thr  following  difference 
equations,  and,  beginning  with  a  particular  set  of  initial  conditions,  advance 
a  problem  In  n  stepe  of  At  over  a  sash  of  I,  Ax  In  accordance  with  sta¬ 
bility  and  accuracy  requirements. 


u“l.u  “-1 
I  I 


a 


»  a  «  n  .  n-i  n-i 

'«♦*  ‘  ?<-i  -«i-i 


m-1  .  n  n+£  .  _ 

1  •  S  *  ui 


V  -  AT 


/  >  V.-) 

Xt+  Ij”  ♦  Xf 


5H" 


U1  * u 


\  <3) 


1- 1 


X 


1-1 


X  ^  ♦  x  n  -  X  n*1 
*1  *1  *1-1 


-  X 
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Tha  main  ft aturs  of  aaeh  problem  It  tbt  outward -moving  shock.  Fig.  1 
compares  tbt  thtek  overpressure  originating  from  tvo  isothermal  sphere  s,  one  Ini¬ 
tially  at  121  atmoephsres,  tha  other  at  2003  atmospheres,  with  a  point-source 

solution  (started  from  an  analytical  form  ralid  for  strong  shocks  and  due  to 

U. 

Taylor  or  von  Neumann.  )  The  isothermal  spheres  were  at  normal  density  origi¬ 
nally  and  eonsefuently  sere  also  at  a  high  temperature.  In  all  cases  the 
finite  source  curve  becomes  Indistinguishable  from  the  point  source  curve  after 
the  shook  from  the  isothermal  sphere  hae  engulfed  a  mate  of  air  ten  times  the 
Initial  mesa. 

A  cold  source  vas  run  In  vhloh  the  sane  Initial  pressure  vas  used  but  the 
teigeraturv  vas  Initially  equal  inside  and  out,  which  mast  a  high  initial 
density.  In  this  case  the  overpressure  did  not  rise  to  the  point-source 
value  until  quite  late  vfaen  -  c  had  fallen  to  shout  a  half  an  atmosphere .  At 
that  point,  however,  the  ease  engulfed  was  again  approximately  tan  times  tbs 

AKD-2806O  "Shock  Hydrodynamics  and  Blast  Waves,”  Los  Alaaos  Scientific 
Laboratory,  28  October  19^. 
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Initial  bmi. 

Fig.  2  shows  some  pressure  prof  lias  as  a  function  of  the  radlua  at 
various  early  times  after  release  of  an  isothermal  sphere  of  normal  density 
gas  at  two  thousand  ataos  Initial  pressure .  Kota  development  of  the  outvard 
shock  and  the  lnvard  rarefaction.  Note  in  particular  the  development  of  an 
inward  shock  folloving  the  rarefaction.  This  shock  does  not  attsdn  a  net 
lnvard  motion  until  the  rarefaction  has  "reflected"  at  the  origin. 

In  Fig.  3  we  view  pres  sure -radius  profiles  near  the  center  and  at  times 
Just  before  and  after  this  lnvard  shock  impinges  on  the  origin.  Arrows  in¬ 
dicate  the  general  movement  with  time  and  hence  the  sequence  far  nearly  con¬ 
stant  time  -intervals. 

81nce  the  shocks  are,  of  necessity,  spread  over  several  mesh  points,  the 
picture  near  a  singular  point  such  as  the  origin  appears  quite  smeared  out 
vith  a  consequent  loss  of  detail. 

Using  mesh  sines  about  one  fourth  the  original  sine,  portions  of  the 
problem  were  re-run.  Admitting  seme  uncertainty  in  positions  and  peak  values 
of  shocks,  a  eoapar  *  r**n  of  the  shock  pressure  for  this  lnvard  moving  shock 
shows  some  appreciable  discrepancy  only  near  the  origin.  After  reflection  the 
discrepancy  decreases  again  in  spite  of  their  quite  different  histories  near 
the  center.  The  dlecrepaney  may  be  attributed  mainly  to  the  difficulty  in 
identifying  the  shook  vhsn  the  rounding  is  eonparabls  to  the  shock  radius. 

This  shock  pressure  is  decreasing  initially  although  moving  lnvard  sines 
the  net  outvard  expansion  behind  the  main  shock  is  dominant.  (Fig.  U). 

If  one  considers  rather  than  the  absolute  pressure,  the  ratio  of  excess 
shock  pressure  over  the  pressure  immediately  in  front  of  the  shock,  a  quantity 
more  characteristic  of  the  strength  of  this  shock,  then  it  le  seen  that  the 
inward  shock  is  increasing  in  strength  but  remains  stationary  until  its  pressure 
ratio  is  greater  than  unity.  (Fig.  3). 

The  reflected  shock  moves  into  a  higher  pressure  and  so  has  a  lover  ratio. 
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The  reflected  shock  dtcrtaMS  about  Inversely  with  the  radial  distance  al¬ 
though  it  hat  a  finite  strength. 

Tbs  density  profiles  at  early  tines  show  a  sharp  compression  spike  at 
the  naln  shock  and  a  rarefaction  eating  immrd  from  the  surface  of  the  initial 
sphere.  (Pig.  6).  The  rare  fast  loo  reduces  the  density  in  ths  sphere  to  a 
very  lov  Talus.  The  inward  shock  is  seen  to  reflect  and  ncrre  outward,  over¬ 
take  the  aaln  shock  and  leave  a  permanent  discontinuity  in  density  (and 
teaperature)  at  ths  point  of  ooalescsnce. 

Ths  temperature  profiles  at  early  tines  (Fig.  7),  show  ths  cooling  effect 
of  ths  rarefaction  and  the  very  modest  heating  done  by  the  aaln  ahock.  Juat 
as  with  ths  density  a  sharp  break  or  discontinuity  is  maintained  at  the  boun¬ 
dary  of  the  initially  heated  sphere.  The  inward  shock  is  seen  in  even  greater 
contrast  here. 

The  dotted  curve  indicates  the  results  with  a  smaller  mssh  else  at  the 
same  time  aa  the  dashed  curve  which  is  for  a  larger  mssh  sits. 

After  sosm  considerable  time  two  similar  problems  with  considerably 
dlffarent  tone  sites  compare  about  as  in  Pig.  8. 

Ths  portion  near  tbs  origin  was  covered  with  a  fine  grid  in  the  case 
illustrated  by  the  dashed  curve.  The  grid  was  coarser  for  the  same  case  near 
ths  main  shock,  however,  which  accounts  for  tbs  poor  reproduction  of  the  second 
shock  at  this  point  high  on  the  back  of  the  main  shock. 

In  addition  to  the  main  shock  and  ths  inward  shock  shown  previously, 
successively  smaller  shocks  develop  from  interaction  with  ths  contact  discon¬ 
tinuity  st  the  edge  of  ths  initial  sphere.  This  boundary  is  marked  by  ticks  in 
Pig.  9*  Pots  ths  reflected  and  the  transmitted  shocks  being  formed. 

The  initial  rarefaction  reflects  at  ths  origin  and  on  moving  outward  in¬ 
teracts  with  this  earns  discontinuity  to  send  beck  a  pulse. 

nils  pulse,  in  turn,  runs  into  ths  outward -moving  second  shook  and  sends 
inward  a  pulse  containing  elements  of  both  rarefaction  and  compression.  In 
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this  ftgurs,  this  pulss  is  shovn  interacting  vlth  ths  origin.  In  tha  uppsr 
curve  s  rarefaction  followed  by  s  shock  is  norlng  inward.  The  middle  curve 
shows  ths  rarefaction  rs fleeted  and  the  shock  still  inward  bound.  Ths  lover 
curve  finds  this  shock  reflected  and  aoving  out. 

Retracing  a  little  in  tins,  the  next  set  of  pressure -prof  lies  (Fig.  10), 
show  ths  aain  shock  decay  with  tine  and  the  noveasnt  of  the  second  shock  after 
reflection  at  tha  center.  In  ths  second  curve  this  second  shock  is  Just 
passing  through  ths  discontinuity,  and  the  minor  pulse  shown  on  ths  previous 
slide  is  about  to  hit  ths  center.  In  ths  third  curve  the  second  shock  is 
riding  up  ths  main  shock  and  a  tertiary  shock  reflected  from  ths  interface  Is 
racing  lnvard.  Ths  fourth  curve  shows  this  third  shock  outward  bound. 

In  ths  first  curve  of  Fig.  II,  tbs  second  shock  is  about  to  close  with 
the  aain  shock  and  ths  third  shosk  is  encountering  ths  interface.  Ths  sec  owl 
curve  shows  the  first  two  shocks  coalesced  and  ths  third  shock  rising  19  ths 
main  shocks. 

Pulses  and  snail  shocks  rattle  back  and  forth  In  thin  hot  interior  sub¬ 
jecting  ths  origin  to  repeated  shooks  and  temperature  rises,  so  that  ths  temp¬ 
erature  gradient  in  the  gas  that  started  out  as  isothsrnal  bee  ones  acre  and 
more  like  that  of  the  point  source  as  tins  goes  on. 
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